Gravitomagnetic effects in Kerr-de Sitter space-time 
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ABSTRACT 

We explicitly worked out the orbital effects induced on the trajectory of a 
test particle by the the weak-field approximation of the Kerr-de Sitter metric. 
It results that the node, the pericentre and the mean anomaly undergo secular 
precessions proportional to k, which is a measure of the non linearity of the theory. 
We used such theoretical predictions and the latest observational determinations 
of the non-standard precessions of the perihelia of the inner planets of the Solar 
System to put a bound on k getting k < 1CT 29 m~ 2 . The node rate of the 
LAGEOS Earth's satellite yields k < 10~ 26 m~ 2 . The periastron precession of 
the double pulsar PSR J0737-3039A/B allows to obtain k < 3 x 10" 21 m" 2 . 
Interpreting k as a cosmological constant A, it turns out that such constraints 
are weaker than those obtained from the Schwarzschild-de Sitter metric. 

Subject headings: Classical general relativity-Approximation methods; equations 
of motion-Experimental tests of gravitational theories-Orbit determination and 
improvement 
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Introduction 



The General Theory of Relativity (GTR) has passed with excellent resu 



observational tests, a s Solar System and binary pulsars observations show (INi 



;s many 



2005 



Will 



2006 



Turyshev 



20081 ). As a matter of fact, the current values of the PPN parameters are in 



agreement with GTR predictions. 

However, some observations seem to question the general relativistic model of 
gravitational interaction on larger scales. On the one hand, the data coming from the 



galactic rotation curves of spiral galaxies (IBinney and Tremaine 



19871 ) cannot be explained 



on the basis of Newtonian gravity or GTR: the existence of dark matter is postulated to 



reconcile the theoretical model with obse rvations; furt 



the mass discrepancy in galactic clusters (j Clowe et al. 



lermore, dark matter can explain 



20061 ) . On the other hand, a lot of 



observations, such as the light cu rves of the type l a 



backg round (CMB) experiments (IRiess et al. 



1998 



supernovas and the cosmic microway e 



Perlmutter et al. 



1999 



Bennet et al. 



20031 ). firmly state that our Universe is now undergoing a phase of accelerated expansion. 
Actually, the present acceleration of the Universe cannot be explained, within GTR, unless 
the existence of a cosmic fluid having exotic properties is postulate, i.e. the so called dark 
energy. 

The cosmological constant is one of the candidat e s to e xplain (in the GTR framework) 



what the dark energy is (see e.g. 



Peebles and Ratra 



(120031 ) and references therein). On 



the other hand, modified gravity models that go beyond GTR have been proposed to 
try to explain current obse rvations and, among these models, f(R) theories of gravity 



(ISotiriou and Faraoni 



20081 ) received much attention in recent years. In these theories the 
gravitational lagrangian depends on an arbitrary function / of the scalar curvature R; 
they are also referred to as "extended theories of gravity", since they naturally generalize, 
on a geometric ground, GTR: namely, when f(R) = R the action reduces to the usual 
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Einstein-Hilbert action, and Einstein's theory is obtained. It is interesting to point out 
that the vacuum solutions of GTR with a cosmological constant are also solutions of f(R) 
gravity vacuum field equations: this is always true in the Palatini formalism, while in metr ic 

am 



1993 



gravity this holds for the solutions wit h constant scalar curvature R ( iFerraris et al. 



Allemandi et al. 



2005; 



Magnano 



19951 ) 



The relevance of the cosmological constant in modern gravitational physics is manifest, 
and it is interesting to focus on the solutions of Einstein's field equations with cosmological 
constant, to investigate its role on different scales. For instance, the Schwarzschild-de Sitter 
metric, which describes a poi nt-like mass in a spa c e-tim e with a cosmologica l cons tant, 



has been recently studie d 



Jetzer and Serend (120061 ) 



by 



K!agramanova et al. 



fl2006h : ISereno and Jetzeri (120061 ) 



Ioriol (l2006bl ). In particular, the Schwarzschild-de Sitter metric 
has been studied to investigate the i nfluenc e of the cosmological con s tant on gravita tional 



lensing in 



Rindler and Ishakl (120071 ): 



Serenol (12008a] ) 



Ruggierol (120071 ) 



Serenol (12008b! ) 



In this paper we are concerned with the Kerr-de Sitter met ric, which descri 



rotat i ng black-hole in a space-time with a cosmolo gical constant (IDemianski 



19731: 



Kerr et al. 



2003 



Kraniotis 



2004 



2005 



Des a 



1973 



Carter 



20071 ). In particular, we want to study the 



gravito-magnetic (GM) effects in Kerr-de Sitter metric. GM effects are due to the rotation 
of the sources of the gravitational field: this gives raise to the presence of off-diagonal terms 
in the metric tensor, which are responsible for a variety of effects concerning orbiting test 
particl es, precessing gyroscopes, moving clocks and atoms and propagating electromagneti c 



waves (IMashhoon et al. 



2001 



Ruggiero and Tartaglia 



2002 



S chafer 



2004 



Mashhoon 



20071 ) 



hey are exp ected in GTR, but are generally very small and, hence, very difficult to detect 



( llorio 



2007aJ). In recent years, th ere have been some attempts to measure the Lense-Thirring 



effect ( iLense and Thirring 



19181 ) with t he LAGEOS and LAGEO S II laser-ranged satellites 



in the gravitational field of the Earth (jCiufolini and Pavlis 



20041 ): the evaluation of the 
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realistic accuracy reached i n such a test and other topics related to it are still matter of 



debate (jlorio 



2006a . 



2007bl ). For other attempts to measure the Lense- 



Thirring eff ect in 



other Solar System scenarios with natural and artificial satelli tes, see (jlorio 



April 2004 the Gravity Probe B spacecraft (jEveritt et al. 



2007ah . In 



20011 ) was launched to ac curate ly 



measure the grav ito-magnetic (and geodetic) precession of an orbiting gyroscope (IPugh 



1959 



Schifj 



I960) in the terrestrial space environment: the final results are going to be 



published. We focus on the GM effects in Kerr-d e Sitte r metric (the GM precession of an 



orbiting gyroscope was investigated by 



Ruggierol (120081 )): in particular we work out the 



GM effects in the weak-field and slow-motion approximation on the orbit of a test particle, 
working out explicitly t he perturbatio ns of the Kepleri an orbital ele ments; furthermore, 



by using the EPM2004 (IPitieva 



2005aT ) and EPM2006 (IPitieval 120081 ) ephemerides, we put 
constraints on the parameter k, which is the cosmological constant in GTR and a measure 
of the non linearity of the theory in f(R) gravity. 



2. Gravito-magnetic field in Kerr-de Sitter metric 



The Kerr-de Sitter metric in the standard Boyer-Lindquist coordinates = (t, p, 6, 



has the form 



ds 2 



2GM p k 
c 2 S ~ 3 



2 ( J_\ h_GMp k 
\Mc) \ c 2 S 3 



J 

Wc 



sin 2 # 



c 2 dt 2 + 



J 

Wc 



sm 2 9cdtd(f) + —dp 2 + —d0 2 - 
A x 



l The space-time metric has signature (1,-1,— 1,-1), greek indices run from to 3, and 
latin ones run from 1 to 3, boldface letters like r refers to three- vectors. 
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(2GMp 



J 

Wc 



sin 



1 + - 



k ( J 
3 \Wc 



P * +[ Wc 



sin 



where 



J 

Wc 

GMp 



cos 2 6 



J 

Wc 



3 \Mc 



(2) 



3 P 



P 2 + 



J 

Wc 



(3) 



The mass of the source is M, while J is its angular momentum (which is perpendicular 
to the 9 = it 1 2 plane); k is the cosmo logical constant in GTR framework, while in f(R) 



theories (jRuggiero 



20081 ) it is a parameter related to the non linearity of the gravity 



lagrangian (namely, when f(R) = R then k — 0). When k — the Kerr-de Sitter metric 
given by eq. ([1]) reduces to the Kerr metric. Other limiting cases can be checked: for 
instance, when J = 0, we obtain the Schwarzschild-de Sitter solution, and when M = J = 
we have the de Sitter space-time. 

In order to study GM effects a weak field approximation eq. (pQ) is sufficient; 
furthermore, it is useful to introduce the isotropic radial coordinate r denotes the defined as 

GM hp 2 



r = p I 1 



c 2 p 



12 



(4) 



where p is the standard Boyer-Lindquist radial coordinate. Then, up to linear terms in 
J^r, kr 2 , the metric is 



GM 



ds 2 = [l- ^ - -r 2 ) c 2 dt 2 - ( 1 



2GM k 



+ 2 



J 

Wc 



c 2 r 3 

2GM kr ( 5GM S 
r 



c 2 r 6 



(dr 2 + r 2 d6 2 + r 2 sin 2 6d(f) 2 ) + 



c 2 r 3 V 2 c 2 / 



sin 2 9d(ficdt. 



(5) 



By differentiating 



y 

— = tan i 



(6) 
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and using 



it turns out that 



x 



COS ( 



sin 



r sin 9 ' 
-y dx + x dy 



(7) 



(8) 



thus, the off-diagonal, i.e. gravito-magnetic, components g^, i = 1,2,3 of the metric 
tensor of the weak-field approximation of the Kerr-de Sitter space-time are, in cartesian 
coordinates, 

J 



9oi 

902 
903 



Mc 
J 
Wc 

0; 



2GM k ( 5GM\ 

3 v + ^ r ; 



y, 



2GM k ( 5GM\ 



c 2 r 3 



(9) 
(10) 
(11) 



In the weak-field and slow-motion linear approximation the spatial c omponents of the 



geodesic equations of motions yielding gravito-magnetic accelerations are (IBrumberg 



19911 ) 



d 2 x l 

~d¥ 



c (djh 0i - dih Qj ) — , i = 1, 2, 3 



(12) 



where h i = g i — r] i = g ot , I = 1, 2, 3. It can be straightforwardly showed that the t erms 



not c ontaining; k yield the usual Lense-Thirring acceleration in cartesian coordinates (jSoffel 



19891 ). The components of the acceleration containing k in cartesian coordinates are 



-4, 



A, 



2GJk 
3c 2 

2GJk 
3c 2 

bGJk 
6c 2 



c 



5 / 

GM + 4 



y 2 + 2z 5 



cr 5 fx 2 + y 2 + 2^ 2 
GM + 4 



z{xy — yx) 



x 



5yz 
4r 3 
5xz 



4r 3 



(13) 
(14) 
(15) 



which can be cast into the vectorial form (IMashhoon 



20071) 



v 

-2- x B, 

c 



(16) 
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where the gravito-magnetic field B is 



„ Jkc ^ 5GJk 

B = ^r-.J + 



j ■ f) f 



3M 12c r 
with J = z; the unit vector in the radial direction is defined as 



(17) 



r 

r = — . 

r 



We notice that the gravito-magnetic field consists of two contributions, the first one that 
is everywhere constant and parallel to J, the second one whose position and directions are 
p osit ion- dep endent . 



Furthermore, by defining the gravito-magnetic potential A. as (jMashhoonl 120071 ) 

9 /c 2 r 5\GJX7- 



A = kr l 



QGM 12 c r 3 



(19) 



and using 



V x (A x B) = (B • V)A - (A • V)B + A(V • B) - B(V • A) 

A x (B x C) = (A • C) B - (A • B) C, 

V x (V>H) = x H + x H, 

V • r = 3, 
_1 r 



(20) 

(21) 
(22) 
(23) 
(24) 



with il> = '. 


/r, it 


( Mashhoon 


2007T 



B = V x A. 



(25) 



In order to calculate the impact of eq. f|T6l) on the orbit of a test particle, let us project 
it onto the radial (f), transverse (t) and normal (n) directions of the co- moving frame 
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picked out by the three unit vectorgj (IMontenbruck and Gill 



2000) 



r = cosw x + cos i sin u y + sin z sin m z, 
t = — sin u x + cos % cos u y + sin i cos m £ 
n = — sini y + cosi £, 



(26) 
(27) 
(28) 



where i is the inclination of the orbital plane to the equator of the central mass and 
u = oj + / is the argument of latitude defined as the sum of the argument of the pericentre 
u, which fixes the position of the pericentre with respect to the line of the nodes, and the 
true anomaly / which reckons the position of the test particle from the pericentre. Thus, 

GJk 



Ay 

An 



A • r 

A-i 

A • h 



6 c 2 



5 + 4 



2 

err 
GM 



0. 



GJk 
~3J~ 



2 

( c r 
''GM 



cos % u, 



sin i sin u u: 



(29) 
(30) 
(31) 



they must be inserted into the right-hand-side of the Gauss equations (IBertotti et al. 
of the variations of the Keplerian orbital elements 

2 



20031) 



da 
~dt 
de 
~dt 
di 
dl 
dfl 
~dt 
du 
~dt 
dM 
~df 



ii 



eA r sin / + A t ( - 



na 
1 



nay/1 — e 2 
1 



A r smf + A t 



A n ( - 1 cos u 



cos/ 



l/_r 
e \ a 



na sin % 



-.A r 



sin u, 



nae 
2 

n A r 

na 



A r cos / + A t ( 1 + - ) sin / 

du 



cost- 



.dn 
~dt' 



dn\ 

, cosi-— 
dt dt J 



(32) 
(33) 
(34) 
(35) 
(36) 
(37) 



2 Here we have chosen the x axis coincident with the line of the nodes, i.e. Q = 0. 
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where a, e, f2 and Ai are the semi-major axis, the eccentricity, the longitude of the ascending 
node and the mean anomaly of the orbit of the test particle, respectively, p = a(l — e 2 ) is 



the semi-latus rectum and n = a/ GMj a? is the un-perturbed Keplerian mean motion. By 
evaluating them onto the un-perturbed Keplerian ellipse 

a(l -e 2 ) 



1 + e cos / 

and averaging^! them over one orbital period Pb of the test particle by means of 

dt (1 - e 2 ) 3 / 2 



(3* 



P b 2tt(1 + ecos/) 2 
it is possible to obtain the secular effects induced by eq. (JTH 



df, (39) 



(a) = 0, (40) 

(e) = 0, (41) 

(i) = 0, (42) 

• \ Jk / 5GM\ 

Q ) = 3m( 1 + ^J' (43) 
Jkcosi / 5GM\ 

<w> = -^( 2+ iaJ- (44) 

/ 3M V c a / 



2^ 



In the calculation we have neglected terms of order 0{e 

The correction APb to the orbital period Pb due to eq. (II 6 1) can be calculated using 
the mean longitude 

X = M + uj + cosi fi. (46) 
For small eccentricities eq. (J35l) -eq. (J37l) yield 

— « n A r - . (47) 



na V a 



3 We used u = f because over one orbital revolution the pericentre u can be assumed 
constant. 
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By using eq. fl29|) it is possible to obtain, for e — >• 

2tt 



ft 



GJA; [Ac 2 a 
1 — — — — I + 5 cos i 



so that 



AP b 



3c 2 na V GM 



2-nJka 2 ( Ac 2 a 



3c 2 M \GM 



5 cosz. 



(48) 



(49) 



We will, now, put constraints on k from the corrections to the standard Newto- 
nian/Einsteinian precessions of the longitudes of the perihelia w = uo + cos 2 Q of the 
inner planets of the Solar System, quoted in Table [fl estimated by E. V. Pit.jeva by fi tting 



more than 400 000 observatio ns of various kinds with the EPM2004 (IPitjeva 



2005a 



and EPM2006 flPitieva 



20081 ) ephemerides. No gravito-magnetic terms of any kind were 
included in the dynamical force models used, so that, in principle, they account for the 
effects investigated by us. For the spin angular momentum of the Sun we will use the valu e 



J Q = (190.0 ± 1.5) x 10 39 kg m 2 s 1 determined from helioseismology (jPijpers 



1998 



2003( 1. 



i.e. independently of the planetary dynamics which we want to test. From 

Jk cos % 



3M 



(50) 



and Table [T] it is possible to obtain 



k < 1 x 10~ 29 m" 2 . 



(51) 



In the case of the laser-ranged LAGEOS satellite (ISmith and Dunn 



19801 ). orbiting 



at about 6000 km above t he Earth's surface, by assuming J ffi = 5.85 x 10 3 3 kg m 2 s 1 



flMcCarthv and Petit 



20041 ) and an uncertainty of the order of 1 cm or less ( ILucchesi 



20071 ) 



in reconstructing its orbit, which translates into an uncertainty in the nodal rate of 5Q ~ 0.1 
milliarcseconds per year, the bound which can be obtained from eq. (1431) is k < 4 x 10 -26 



m 
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Concerning the double pulsar system PSR J0737-3039A/B (IBurgav et all 120031). by 



assum ing for the moment of inertia of A t he value I 



20051 ). since its rotational period is 22 ms (IKramer et al. 



0f_kg m (lLorimer and Kramer 



20061 ) its angular momentum can 



be evaluated as Ja = 2.8 x 10 40 kg m 2 s 1 . The overall uncertainty (including also the 
mism odelling in the usual 1PN term) in the periastron precessi on amounts to 0.03 deg 



yr 1 (jlorio 



20091 ). so that eq. (jSJ) and the system's parameters (IKramer et al.l 120061 ) yield 



k < 3 x 10" 21 m" 2 . 



In order to evaluate such bounds it is useful remember that studying the Schwar z schild 



de Sitter non-gravitomagnetic 


precession various authors ( 


Kagramanova et al. 


20 


06; 


Iorio 


2006b 




Sereno and Jetzer 


2006 




Jetzer and Sereno 


2006 




Ruggiero and Iorio 


2007 


) have 



obtained bounds on the cosmological constant that range from 10~ 40 m -2 to about 10~ 42 
m -2 ; furthermore, the data in the general relativistic ACDM (cosmological constant plus 
cold dark matter) model suggest that a value of the cosmological const ant of 10~ 52 m~ 2 is 



1996|). 



required in order to explain the accelerated expansion of the Universe ( IKochanek! 
As a consequence, the bounds that we have just obtained are not competitive. The same 
conclusions apply to Palatini f(R) gravity: in fact when the formation of large-scale 
structure is studied in this context (i.e. when the non linear part of f(R) is supposed to 
drive cosmic ac celeration ) the r esults are practically indistinguishable from the ACDM 



model; see e.g. iKoivistol (120061 ) who used large scale structure cosmological data from 
Sloan Digital Sky Survey (SDSS). However, these tight constraints on f(R) gravity are 
loosened when dark m atter with inherent stresses (generalized dark matter, GDM) is 



allowed (IKoivisto 



20071). 



3. Discussion and Conclusions 



We explicitly worked out the effects induced on the orbit of a test particle by the 
weak-field approximation of the Kerr-de Sitter metric, which is a solution of the vacuum 
field equations both in GTR and in f(R) gravity. It turns out that the semi-major axis, 
the eccentricity and the inclination do not experience secular, i.e. averaged over one orbital 
period, changes; instead, the longitude of the ascending node, the argument of pericentre 
and the mean anomaly undergo secular precessions. Interestingly, all such effects consist 
of two kinds of contributions. The first type is given by terms proportional to GMc~ 2 , 
which vanish in the limits c — > oo, G — > 0, M — > 0. Instead, the second kind consists of 
terms proportional to Jk/M, which are independent of the speed of light c, the constant of 
gravitation G and the source's mass M, so that they do not vanish in the limits for c — ■> oo, 
G — > and M — > 0. Concerning the dependence on the orbital geometry, both kinds of 
effects depend on the inclination and vanish for polar orbits; while the 0(c~ 2 ) terms depend 
also on the size of the orbit through the semi-major axis, it is not so for the O(JkfM) 
ones which are, indeed, independent of it. Then, we compared our predictions to the latest 
observational determinations of the corrections to the standard Newtonian/Einsteinian 
precessions of the perihelia of the inner planets of the Solar System obtaining the constrain 
k < 10" 29 m~ 2 . The node of the terrestrial LAGEOS satellite yields k < 10~ 26 m" 2 , 
while the bound from the periastron of the double pulsar system PSR J0737-3039A/B is 
k < 3 x 10~ 21 m~ 2 . Such bounds are not competitive with the ones which can be obtained 
from the Schwarzschild-de Sitter non-gravitomagnetic precessions (k < 9 x 10~ 43 m -2 from 
Solar System data) and from those deriving from cosmological observations. 



-14- 



Table 1: Inner planets. First row: estimated peri helion extra-pr ecessions in 10 



-4 // 



cy" 



arcseconds per century), from Table 3 of (IPitjeva 



cy 



-i /// 



2005bl ) (apart from Venus). The 



quoted errors, in 10 4 " cy , are not the formal ones but are realistic. The formal errors are 
quoted in square brackets (E.V. Pitjeva, personal communication to L.I., November 2005). 



The units are 10 4 " cy l . Second row: semi-ma 



formal errors are in Table IV of (IPitjeva 



or axes, in Astronomical Units (AU). Their 



2005al ). in m. Third row: eccentricities. Fourth 



row: orbital periods in years. The result for Venus have been recently obtained by including 
the Magellan radiometric data (E.V. Pitjeva, personal communication to L.I., June 2008). 



Mercury 



Venus 



Earth 



Mars 



(Aw) (KT 4 " cy^ 1 ; 
a (AU) 

e 

P (yr) 



-36±50[42] -4±5[1] -2±4[1] 1 ± 5[1] 

0.387 0.723 1.000 1.523 

0.2056 0.0067 0.0167 0.0934 

0.24 0.61 1.00 1.88 
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